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Abstract 

In many problems of classical analysis extremal configurations appear 
to exhibit complicated fractal structure. This makes it much harder to 
describe extremals and to attack such problems. Many of these problems 
are related to the multifractal analysis of harmonic measure. 

We argue that, searching for extremals in such problems, one should 
work with random fractals rather than deterministic ones. We introduce a 
new class of fractals random conformal snowflakes and investigate its prop- 
erties developing tools to estimate spectra and showing that extremals can 
be found in this class. As an application we significantly improve known 
estimates from below on the extremal behaviour of harmonic measure, 
showing how to constuct a rather simple snowflake, which has a spectrum 
quite close to the conjectured extremal value. 



1 Introduction 

It became apparent during the last decade that extremal configurations in many 
important problems in classical complex analysis exhibit complicated fractal 
structure. This makes such problems more difficult to approach than similar 
ones, where extremal objects are smooth. 

As an example one can consider the coefficient problem for univalent func- 
tions. Bieberbach formulated his famous conjecture arguing that the Koebe 
function, which maps a unit disc to a plane with a stright slit, is extremal. The 
Bieberbach conjecture was ultimately proved by de Branges in 1985 [3J, while 
the sharp growth asymptotics was obtained by Littlewood [B] in 1925 by a much 
easier argument. 

However, coefficient growth problem for bounded functions remains widely 
open, largely due to the fact that the extremals must be of fractal nature (cf [2]). 
This relates (see [T]) to a more general question of finding the universal multi- 
fractal spectrum of harmonic measure defined below, which includes many other 
problems, in particular conjectures of Brennan, Carleson and Jones, Kraetzer, 
Szego, and Littlewood. 

In this paper we report on our search for extremal fractals. We argue that 
one should study random fractals instead of deterministic ones. We introduce 
a new class of random fractals, random conformal snowflakes, investigate its 
properties, and as a consequence significantly improve known estimates from 
below for the multifractal spectra of harmonic measure. 
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1.1 Multifractal analysis of harmonic measure 

It became clear recently that appropriate language for many problems in geo- 
metric function theory is given by the multifractal analysis of harmonic measure. 
The concept of multifractal spectrum of a measure was introduced by Mandel- 
brot in 1971 in [51 [TU] in two papers devoted to the distribution of energy in 
a turbulent flow. We use the definitions that appeared in 1986 in a seminal 
physics paper [1] by Halsey, Jensen, Kadanoff, Procaccia, Shraiman who tried 
to understand and describe scaling laws of physical measures on different fractals 
of physical nature (strange attractors, stochastic fractals like DLA, etc.). 

There are various notions of spectra and several ways to make a rigorous 
definition. Two standard spectra are packing and dimension spectra. The pack- 
ing spectrum of harmonic measure u> in a domain Q with a compact boundary 
is defined as 



7r n (t) = sup lq : Vtf > 3 S - packing {B} with diam(S)*w(S) <? > 1 



where (5-packing is a collection of disjoint open sets whose diameters do not 
exceed S. 

The dimension spectrum which is defined in terms of harmonic measure to 
on the boundary of Q (in the case of simply connected domain £1 harmonic 
measure is the image under the Riemann map <p of the normalised length on 
the unit circle). Dimension spectrum gives the dimension of the set of points, 
where harmonic measure satisfies a certain power law: 



Here dim stands for the Hausdorff or Minkowski dimension, leading to possibly 
different spectra. The restriction a > 1/2 is due to Beur ling's inequality. Of 
course in general there will be many points where measure behaves differently 
at different scales, so one has to add lim sup's and liminf's to the definition 
above - consult [7] for details. 

In our context it is more suitable to work with a modification of the packing 
spectrum which is specific for the harmonic measure on a two dimensional simply 
connected domain f2. In this case we can define the integral means spectrum as 



where is a Riemann map from the complement of the unit disc onto a simply 
connected domain SI. 

Connections between all these spectra for particular domains are not that 
simple, but the universal spectra 





P<ji{t) '.= lim sup 

r- »1+ 



logj^l^re*)!^ 
|log(r-l)| 



, teK, 



U(t) = supn(t), F(a) = sup /(a), and B(t) = svp0(i) 




2 



are related by Legendre-type transforms: 



F(a) = inf (allft) +t), a> 1 , 

0<t<2 

H(t) = supf^M^lV 0<2<2, 
U(t) = B(t) -t + 1 . 
See Makarov's survey [7] for details. 

1.2 Random fractals 

One of the main problems in the computation of the integral means spectrum 
(or other multifractal spectra) is the fact that the derivative of a Riemann map 
for a fractal domain depends on the argument in a very non regular way: <j)' is a 
"fractal" object in itself. We propose to study random fractals to overcome this 
problem. For a random function (j> it is natural to consider the average integral 
means spectrum: 

p(t) = supj/3 : j{r- l f^J a E [|/'(re*)|*] d0dr = oo J 

= inf j/3 : J(r- 1 ) /3_1 ^ E ^/'(re* 9 )!'] dOdr < oo 

The average spectrum does not have to be related to the spectra of a particular 
realization. We want to point out that even if <f> has the same spectrum a.s. it 
does not guarantee that f3(t) equal to the a.s. value of (3(t). Moreover, it can 
happen that (3 is not a spectrum of any particular domain. 

But one can see that (3{t) is bounded by the universal spectrum B(t). Indeed, 
suppose that there is a random / with (3(t) > B + e, hence for any r there 
are particular realizations of / with J\f'(z)\d9 > (r — l)~ B ~ e / 2 . Then by 
Makarov's fractal approximation _7j there is a (deterministic) function F such 
that /3f(*) > B(t) which is impossible by the definition of B(t). 

For many classes of random fractals E|0'| 4 (or its growth rate) does not 
depend on the argument. This allows us to drop the integration with respect to 
the argument and study the growth rate along any particular radius. Perhaps 
more importantly E\<j)'\ is no longer a "fractal" function. 

One can think that this is not a big advantage compared to the usual integral 
means spectrum: instead of averaging over different arguments we average over 
different realizations of a fractal. But most fractals are results of some kind 
of an iterative construction, which means that they are invariant under some 
(random) transformation. Thus E|0'|* is a solution of some kind of equation. 
Solving this equation (or estimating its solutions) we can find (3(t). 

In this paper we want to show how one can employ these ideas. In the 
Section [2] we introduce a new class of random fractals that we call random con- 
formal snowflakes. In the Section [3] we show that (3(t) for this class is related to 
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the main eigenvalue of a particular integral operator. We also prove the fractal 
approximation for this class in the Section |4j In the Appendix [5] we give an 
example of a snowflake and prove that for this snowflake (3(1) > 0.23. This sig- 
nificantly improves previously known estimate B(l) > 0.17 due to Pommerenke 
PI- 

2 Conformal snowflake 

The construction of our conformal snowflake is similar to the construction in 
Pommerenke's paper [11] . The main difference is the introduction of the ran- 
domness. 

By £' we denote a class of all univalent functions <p : B_ — > B_ such that 
4>(oo) = oo and 4>'(oo) G M. Let <f> G £' be a function with expansion at infinity 
<p(z) = b\Z + . . . , by cap (f> = capQ we denote the logarithmic capacity of </> 
which is equal to log We will also use the so called Koebe n-root transform 
which is defined as 

(K<f>)(z) = (K n cf>)(z) = ^/H^). 

It is a well known fact that the Koebe transform is well defined and K<f> € 
It is easy to check that Koebe transform divides capacity by n and that the 
capacity of a composition is the sum of capacities. 

First we define the deterministic snowflake. To construct a snowflake we 
need a building block and an integer k > 2. Our snowflake will be the 

result of the following iterative procedure: we start with the building block and 
at n-th step we take a composition of our function and the fc™-root transform 
of the rotated building block. 

Notation 1. Let <f> € £' and 8 £ [0, 2w]. By <j>e(z) we denote the map whose 
range is the rotation of that for cf>, namely e* 6 '</>(ze -1 ' ) . 

Definition 1. Let <fi 6 E', k > 2 be an integer number, and {9 n } be a sequence 
of numbers from T. Let fo(z) = 4>e {z) and 

/»(*) = fn-i(K k ^ en (z)) = 9o (^ fc (. . . ^(z^) ...). 

The conformal snowflake / is the limit of /„. For simplicity S = C \ /(©_) and 
g = f^ 1 are also called a snowflake. 

Sometimes it is easier to work with a slightly different symmetric snowflake 

/„(*) = ti{ k (. . . < 7 > fe ") • • •) = $1 o • • • o $„(z), 

where $j = if fe j c/dg . . There are two equivalent ways to construct the symmetric 
snowflake from the usual one. One is to take the Koebe transform Kf n , another 
is to start with fo(z) — z. It is easy to see that /„ = $o 0, "° 

How this snowflake grows? This is easy to analyse looking at the evolution 
of f n . At every step we add k n equidistributed (according to the harmonic 



4 



Figure 1: Images of /i, f 2 , h, and / 4 



measure) small copies of the building block. But they are not exact copies, they 
are distorted a little bit by a conformal mapping. 

Figures [T] and [2] show images of the first four functions / and / with k = 2 
and the building block is a slit map (which adds a straight slit of length 4). 

Lemma 2.1. Let f n = cf>g (f n (z)) be the n-th approximation to the snowflake 
with a building block <f> and k > 2. Then cap (/„) and cap (/„) are bounded by 
(and converge to) cap ((f>)k/(k — 1) and cap ((/>)/ (k — 1). 

Proof. This lemma follows immediately from the standard facts that 

cap (/ o g) = cap (/) + cap (g) , 
cap (K n f) = cap (f)/n. 

□ 

Theorem 2.2. The conformal snowflake is well defined, namely let f n be an 
n-th approximation to a snowflake with a building block (f> and k > 2. Then 
there is f € £' such that f n converge to f uniformly on every compact subset 
of®- 

Proof. Fix e > 0. It is enough to prove that /„ converge uniformly on D e = 
{z : \z\ > 1 + e}. Suppose that m > n so we can write f m — f n o $> ntm where 
$n,m = ®n+l o • • • o $ m and 

\fn(z) ~ f m (z)\ = \fn(z) ~ fn($n,m(z))\ < max |/^(C)| \z - $ n ,m(z)\ . 
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Figure 4: The image of a small boundary arc under with three Green's lines. 

By the Lemma 12. II cap (/„) is uniformly bounded, hence by the growth theorem 
the derivative of /„ is uniformly bounded in B e . Thus it is enough to prove that 
$n,m(z) converge uniformly to z. 

Let 4>{z) = b\z + . .. at infinity, then <f> n (z) — b\^ k z + . . . Functions $„. m 
have expansion 

b k 1 ~ n+ - +k ~ m z + --- = b i ; i ' m) z + ... 

Obviously, b^" — > l as n -> oo. This proves that ^ n<m (z) — > z uniformly on 
D e hence /„ converge uniformly. Uniform limit of a functions from £' can be 
either a constant or a function from Since cap (/„) is uniformly bounded 
the limit can not be a constant. □ 

Definition 2. Let (f> <E £' and k > 2 be an integer number. The random 
conformal snowflake is a conformal snowflake defined by 0, k, and {# ra }, where 
n are independent random variables uniformly distributed on T. 

Theorem 2.3. Let <p £ S', k > 2 be an integer number, and ip — </> -1 . Let f 
be a corresponding random snowflake and g = f . Then the distribution of f 
is invariant under the transformation S' x T h> S' which is defined by 

(f,6)^MK k f). 
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In other words 

f(z) = M(K k f)(z)) = Mf 1/k (z k )), (1) 
g(z) = (K k g)(Mz))=g 1/k (^(z)), (2) 

where 9 is uniformly distributed on T. Both equalities should be understood in 
the sense of distributions, i.e. distributions of both parts are the same. 

Proof. Let / be a snowflake generated by {#„}. The probability distribution 
of the family of snowflakes is the infinite product of (normalised) Lebesgue 
measures on T. By the definition 

/(*)= lim eo (^ fc (...^ fc (z fc ")...) 

and 

M(K k f)(z))= lim e (^ o /fc (...^/ fc (z fc " +1 )...), 

hence <j>e {Kkf) is just the snowflake defined by the sequence 9, 9q, 6\, . . . . So the 
transformation f(z) i— > <pe({Kkf){z)) is just a shift on the [— tt, 7r] N . Obviously 
the product measure is invariant under the shift. This proves stationarity of /. 
Stationarity of g follows immediately from stationarity of /. □ 

There is another way to think about random snowflakes. Let M. be a space 
of probability measures on S'. And let T be a random transformation / 
4>e(Kkf), where 9 is uniformly distributed on [tt,tt]. Obviously T acts on M. 
The distribution of a random snowflake is the only measure which is invariant 
under T. In some sense the random snowflake is an analog of a Julia set: it 
semi-conjugates z k and ipg. 




The random conformal snowflakes are also rotationally invariant, the exact 
meaning is given by the following theorem. 
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Theorem 2.4. Let <fi G X', k > 2 and <? 6e </ie corresponding snowflake. Then 
g is rotationally invariant, namely g(z) and e %UJ g(eT lu z) have the same distri- 
bution for any u>. 

Proof. Let g n {z) be the nth approximation to the snowflake defined by the 
sequence of rotations 9o, . . . ,6 n . We claim that g(z) = e %0J g{e~ %0J z) is the ap- 
proximation to the snowflake defined by 80, . . . ,9 n where 6j = 9j +ioki (we add 
arguments mod 27r). 

We prove this by induction. Obviously this is true for go. Suppose that it 
is true for g n -i- By the definition of g n and assumption that g n -\{e~ w z) — 
e~ luJ g n -\{z) we have that 

e w g n {e' lu} z) = e^e*"/*" V> 1/fe " (e _ *"ff£" i(e _<a, «)) 
= e »„/fc-^i/fc» (e w„^» = ^ (z) . 

Obviously n are also independent and uniformly distributed on T], hence 
g n has the same distribution as g n . □ 

Corollary 2.5. The distributions of \g(z)\ and \g'(z)\ depend on \z\ only. The 
same is true for f. 

3 Spectrum of a conformal snowflake 

As we discussed above, for random fractals it is more natural to consider the 
average spectrum 0{t) instead of the usual spectrum 0(t). We will work with 
0(t) only and "spectrum" will always mean 0{t). 

Notation 2. We will write L for the class of functions on (1, 00) that are bounded 
on compact sets and integrablc in the neighbourhood of 1. In particular, these 
functions belong to L 1 ^, R] for any 1 < R < 00. 

Let F(z) = F(\z\) = F{r) = E [ \g'{r)/g{r)\ T \og a \g{r)\] where r = 2 - t and 

(7 = 0-1. 

Lemma 3.1. The 0(t) spectrum of the snowflake is equal to 

mf{0 : F(r) e L} . 
Proof. By the definition is the minimal value of such that 



{r-lf-^fire^Wdedr 

is finite. We change variable toto = f(z) — fire" 10 ) 

E[|/'(z)| t ](|z|-l)' 3 - 1 dm(z) 



J J E [|/'(re i9 )|*] (r - lf^drdO = J 

-I 



E 



r 

i-t(\„(„,,\\ -\\0-~t 



W(w)\ 2 -\\g(w) \ - 1) 

\gM\ 



dm , 
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where m is the Lebesgue measure. Note that \g\ is uniformly bounded and g is 
rotationally invariant, hence the last integral is finite if and only if 



\g'(r)\mg(r)\-iydr<™. 



Since 1 < \g\ is uniformly bounded we have that |<7 / (t')| t (|<?(r)| — l) a is compa- 
rable up to an absolute constant to 



\g'(r)\y 

\9(r)\ 



log CT |ffWI- 



Lemma 3.2. If F G L then it is a solution of the following equation: 

d0_ 
2?r' 



F(r) = - J F(|^V <9 )IW*~Ve < V(re <9 )r; 



□ 



(3) 



Proof. By the Theorem 12.31 a(z) and g 1 ^ k {4'g{ z )) have the same distribution, 
hence 



F{r) 



E [\g\r)/g{r)\ T log* \g[r)\] 

(s^WW))' 



E 



E 



s 1/fe WW) 



»(#)) 



log CT |5 1/fc (^(r))| 



log CT ( 5 (^(r))) 



where has a uniform distribution. The expectation is the integral with respect 
to the joint distribution of g and 9, since they are independent this joint dis- 
tribution is just a product measure. So we can write it as the double integral: 
first we take the expectation with respect to the distribution of g and than with 
respect to the (uniform) distribution of 6 



F(r) = 



log"(ff(i#(r))) 



* s 



lo^ (g(^(r)))dg 



2tt 



The inner integral is equal to F{tK{r)) = F(ip (e l6 r)) by the definition of F, 
hence 



F(r) 



— j FW>*(e"r)W(.e»r)il> k -He»r)\ r — 



which completes the proof. 



□ 
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This equation is the key ingredient in our calculations. One thinks about 
F as the main eigenfunction of an integral operator. Hence the problem of 
finding the spectrum of the snowfiake boils down to the question about the 
main eigenvalue of a particular integral operator. Usually it is not very difficult 
to estimate the latter. 

This justifies the definition: 

/7T jn 
f(\iP k (re i8 )\) l^-WCre")!^ ■ (4) 

Using this notation we can rewrite ([3]) as 

k^F = QF. 

Note that this is in fact an ordinary kernel operator, \ip\ is a smooth function 
of 6, hence we can change the variable and write it as an integral operator. As 
mentioned above, the study of a F is closely related to the study of operator Q 
and its eigenvalues. And our estimate of the spectrum is in fact estimate of the 
main eigenvalue. 

3.1 Adjoint operator 

First of all we want to find a formally adjoint operator. Let v be a bounded 
function and R > 1 such that Dr C tp k (Dfj) where Dr — {z : 1 < \z\ < R}. 



Qf{r)v{r)dr = "( r ) k J q /(^ fc (re i9 )) |^'(re ie )^ 1 (re ie )r^dr 
ID R \z\ 2tt 

where dm is the Lebesgue measure. Changing a variable to w = ip k (z) we get 



/ 



v{(/>(w 1/k )) 1 , 1/k i/fc_i |2 _Tj / 



f* ^ k rv{ct){r 1 / k e ie l k )) il^K^ 9/fc 2 _ r _dg_ 



1 Jo 



So we define another operator 

f 2x ,Arh(r 1 / k p i9 



Changing 2 — r to t we can rewrite (J5)) as 
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The inequality above can be written as 

Qf{r)u(r)dr > / f{r)Pu{r)dr. (7) 



We would like to note that for R = oo there is an equality since ip K (H)_) covers 
ID>_ exactly k times. In this case 

/oo 
f{r)Pv{r)dr. (8) 

so operators P and Q are formally adjoint on [l,oo). 

Lemma 3.3. Operator Q = Q{t) acts on L if J \(f>'(re te )\d8 is bounded. 1ft > 1 
then it also acts on L (1, oo). 

Proof. Let v = 1 in ©. Then 

/w£ 1 ^ I i«^)i , S- 

Let r < R, then 

i (fc-i)« 
r k 1 (fe— t)t 
< R k 

|0( r i/fc e ie)| v 

so the second integral is bounded since / |^>'|*d0 is bounded. This proves that 
Qf is in L. To prove that it acts on oo) we should consider the large values 
of r. At infinity </>(z) = cz + . . . and </>'(.z) = c + . . . , hence 

. _ (k—iyt n (fc-i)t , . 



/// 1/fc iflMt r fc c i i 

W r L l K p w \ r Brl !_!_ — rrmst r 1 fc fc 



'(r v e")l ~ = const r 1 

= const r" 



|^( r l/fe e ^)| lr v " IcIrVfc 

■=-i)(i-t) 



Thus the second integral is comparable (up to a universal constant) to r~ ^ , 
so it is bounded if and only if t > 1. □ 

Remark 1. Note that the assumption on the integral of |0'| is just a bit stronger 
than (3<j,(t) — 0. We restrict ourselves to the building blocks that are smooth 
up to the boundary, for such building blocks this assumption is always true. 
Condition t > 1 is technical and due to the behavior at infinity which should be 
irrelevant. Introducing the weight at infinity we can get rid of this assumption. 
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Next we want to discuss how eigenvalues of P and Q are related to the 
spectrum of the snowfiake. If F is integrable then it is a solution of ([3]) and 
using (O we can write 

*m»(r) = f" > j\(rMr)-^L. (9) 

Suppose that t is fixed. Let us fix a positive test function v. If Pv(r) > v(r)k a+1 
then we arrive at contradiction, this means that F(r) for this particular pair of 
t and a can not be integrable. Using this fact we can estimate j3(t) from below. 
Hence any positive v gives the lower bound of the spectrum. 

f3(t)> min log ( ] /log A;. (10) 



l<r<R \ v(r) 

Obviously, the best choice of v is an eigenfunction of P corresponding to the 
maximal eigenvalue. This proves the following lemma: 

Lemma 3.4. Let A be the maximal eigenvalue of P (on any interval [1, R] such 
that D R C ip k (D R )) then (3(t) > log A/ log k. 

4 Fractal approximation 

In this section we prove the fractal approximation by conformal snowflakes. 
Namely we show that for any t one can construct a snowfiake with building 
block which is smooth up to the boundary and (3{t) arbitrary close to B(t). The 
proof of this theorem is similar to the proof of the fractal approximation for 
standard snowflakes but it is less technical. 

Theorem 4.1. For any e and t there are a building block <f> G E'cap C°°({|z| > 
1}) and a positive integer k that define the snowfiake with (3(t) > B(t) — e. 

We will use the following lemma. 

Lemma 4.2. For any e > 0, t 6 M there is A > such that for any 5 > there 
is a function (f> € S'cap C°° such that 

\<j)'(re te )\ d9> 4 - 

for S > r — 1 . Moreover, capacity of <f> is bounded by a universal constant that 
does not depend on 6. 

Proof. There is a function / with f3f(t) > B(t) — e. Hence there is a constant 
A such that 

" \f , (re l9 )\ t d0> A(r-l)- B ^+ 2 *. 

The only problem is that this function is not smooth up to the boundary. Set 
4>(z) — f(sz). Obviously, <p(z) =i f(z) as s — > 1. If we fix a scale 5 then there is 
s sufficiently close to 1 such that J \<j>' /(j^dO > A5~ B{t)+2li /2. But for r < 1 + 5 
the integral can not be smaller by the subharmonicity. □ 
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Proof of the Theorem \4-l\ It is easy to see that 

cap (/„) < cap (/) = cap (<(>)/ (1 - 1/k) < leap (0), 

hence cap (/„) and |/ n (z)| for \z\ < 2 are bounded by the universal constants 
that depend on capacity of <p only and do not depend on k. It also follows that 
\Kkfn{z)\ < 1 + c/k for \z\ < 2 and c depending on cap {<f>) only. 

Let us fix t and let <f> be a function from the Lemma [4.21 for 6 = c/k. By 
/(/, S) we denote 

f Ifire^fdO, 

J —IT 

where r = exp(£). 

The k-TOot transform changes integral means in a simple way: 

t 

r^-^aie. 



fit k ikB\ 

I(K k f,5/k)= ' 



j(k-l)/k^ r k e ik6^ 

As we mentioned before, the capacity of the snowflake is bounded by the uni- 
versal constant, hence |/| can be bounded by a universal constant. Thus 

I(K k f,S/k)>constI(f,S). 

The function f n +i is a composition of a (random) function tpg with K f n . 
The expectation of I(/„+i, l/k n+1 ) conditioned on /„ is 



E [7(/ n+1 , l/k n+1 ) | /„] = / j \^ g {K k f n {re^))\%K k f n )\re^)\ t d£,d9 

|(K fc /„)'(re l «)|* / Wie-^KkUre^dOdli, 



where r = exp(l/fc n+1 ). We know that \Kkf n (re )\ — 1 < c /k- By our choice 
of 6 



\<fi'(\K k f n (re i t)\e- i6 )\ t d6>A 



k \ B{t)-e 



So 



k \ B(t)-e 



E [I(f n+ i, l/k n+1 )] > A(-) E[l(Kf n ,l/k n+1 )} 



h \ B(t)-e 

> A[-j const E [/(/„, l/k n )] 



k \ "(B(i)-e) 



Applying this inequality ri times we obtain 

E [/(/„, l/k n )} > const" 
So 

logE[I(/ !/*«)] > 
n log fc 

for sufficiently large fc. This completes the proof. □ 
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5 Appendix: example of an estimate 



The main purpose of this section is to show that using conformal snowflakes it 
is not very difficult to find good estimates. Particularly it means that if one of 
the famous conjectures mentioned in the introduction is wrong, then it should 
be possible to find a counterexample. 

In this section we will give an example of a simple snowflake and estimate 
its spectrum at t = 1 . We could do essentially the same computations for other 
values of t, but B(l) is of special interest because it is related to the coefficient 
problem and Littlewood conjecture (see pQ for details). 

As a building block we use a very simple function: a straight slit map. We 
use the following scheme: first we define a building block and this gives us 
the operator P. By (|10p any positive function v gives us an estimate on the 
spectrum. To choose we find the first eigenvector of discretized operator 
P and approximate it by a rational function. We compute Pv using Euler's 
quadrature formula and estimate the error term. The minimum of Pv jv gives 
us the desired estimate of j3(l). For t = 1 we give the rigorous estimate of 
the error term in the computation of Pv/nu, for other values of t we give 
approximate values (computed with less precision) without any estimates of the 
error terms. 



5.1 Single slit domain 

We use a straight slit functions. First we define the basic slit function 



*(*, o = m = „ [-^ /TTWmTW - ) . (id 

where s is a constant close to 1, [i\ and /X2 are the Mobius transformation that 
maps B_ onto the right half plane and its inverse: 

z-1 
A^) = — 

M*) = — 

We also need the inverse function 

1>{z,l)=1n{z) = <l>(z,i)- 1 - (12) 

The function <j> first maps D_ onto the right half-plane, than we cut off a straight 
horizontal slit starting at the origin and map it back. The image <^(B_) is B_ 
with a horizontal slit starting from 1. The length of the slit is The derivative 
of a slit map has a singularities at points that are mapped to 1. But if we take 
s > 1 then these singularities are not in D_. We set s = 1.002. 

We study the snowflake generated by <fi(z) = ^73(0) with k — 13 (numbers 
13 and 73 are found experimentally). Figure [3] show the image of the unit circle 
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under fa . The Figure 2] shows the image of a small arc under and three 
Green's lines. 

First we have to find the critical radius R such that Dr C ip (Dr). By 
symmetry of </>, the critical radius is the only positive solution of 

^ k {x) = x. 

This equation can not be solved explicitly, but we can solve it numerically (we 
don't care about error term since we can take any greater value of R). The 
approximate value of R is 76.1568. To be on the safe side we fix R = 76.2. 
The disc takes just a small portion of i/j(Dr) which means that there is a huge 
overkill in the inequality ([7]). 

By (jlOp any positive function v gives a lower bound of spectrum. And this 
estimate is sharp when v is the main eigenfunction of P. So we have to find an 
"almost" eigenfunction of P. 



5.2 Almost eigenfunction of operator P 

Even for such a simple building block we can not find the eigenfunction explicitly. 
Instead we look for some sort of approximation. The first idea is to substitute 
integral operator P by its discretized version. Here we use the most simple and 
quite crude approximation. 

Choose sufficiently large N and M. Let r n = 1 + (R — l)n/N and 6 m — 
2wm/M. Instead of P we have an N x N matrix with elements 

^ ~V n \<t>(rl' k e^)\M' 

where summation is over all indexes m such that r n i is the nearest point to 
\<j>(rn e %Bm )\. This defines the discretized operator P/y Let An and Vn be 
the main eigenvalue and the corresponding eigenvector. A priori, X N should 
converge to k^*', but it is not easy to prove and not clear how to find the rate 
of convergence. But this crude estimate gives us the fast test whether the pair 
<j) and k defines a snowflake with large spectrum or not (this is the way how we 
found k = 13 and I = 73). 

Instead of proving convergence of Ajv and estimating the error term we will 
study Vjv which is the discrete version of the eigenfunction. We approximate 
V/v by a rational function of a relatively small degree (in our case 5), or by any 
other simple function. In our case we find the rational function by the linear 
least square fitting. In any way we get a nice and simple function v which is 
supposed to be close to the eigenfunction of P. 

We would like to note that procedure, by which we obtained v, is highly 
non rigorous, but that does not matter since as soon as we have some explicit 
function v we can plug it into P and get the rigorous estimate of [3. 

In our case we take N = 1000 and M = 500. The logarithm of the first 
eigenvalue is 0.2321 (it is 0.23492 if we take s = 1). Figure [5] shows a plot with 
coordinates of the first eigenvector. 
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Figure 5: Coordinates of the eigenvector 

We scale this data from [1, 1000] to the interval [1, R] and approximate by a 
rational function v. 

v(x) = (7.1479 + 8.9280a; - 0.07765a: 2 + 1.733 x 10" 3 x 3 - 
2.0598 x 10~ 5 .x 4 + 9.5353 x 10~V)/(2.7154 + 13.2845x). 

Figure [6] shows the plot of v. 

1 

. 9 
. 8 
0.7 
. 6 

6 10 20 30 40 50 60 70 
Figure 6: An "almost" eigenfunction 





5.3 Estimates of derivatives 

To estimate /3(1) we have to integrate i/(|<^|)|</>'|/|(/>|. It is easy to see that the 
main contribution to the derivative is given by a factor |0'|/|0|. Assume for a 
while that s = 1. The fraction \<p' (z) / (f>(z)\ can be written as 

(13) 



\z\,J\(z - zi)(z - z 2 ) 
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where 

^ = -5033-gWg ^ _ Q 894756 _ Q 
5625 



z 2 = - 503 lt o 2 f WU * -0.894756 + 0.4465561 
5625 

Singular points zi and z^ are mapped to 1 and 4>' has a square root type 
singularity at these points. They will play essential role in all further calcula- 
tions. We introduce notation Z\ = x + iy and Z2 = x — iy, for z we will use polar 
coordinates z = re 1 ® . 

We compute the integral of / = \v((f))(f)' / (f>\ using the Euler quadrature for- 
mula based on the trapezoid quadrature formula 

- ■ n-l 

f(x)dx a S?(f) = S e (f) -5>2fee 2fc (/ (2fe - 1) (2 7 r) - /^(O)) , 

fc=i 

where S e (f) is a trapezoid quadrature formula with step e and jk — B k /k\ were 
Bk is the Bernoulli number. The error term in the Euler formula is 



72„ max 



/( 2 ™) e 2 "27r. (14) 



In our case function / is periodic and terms with higher derivatives vanish. This 
means that we can use (|14p for any n as an estimate of the error in the trapezoid 
quadrature formula. 

Function <f> has two singular points: z\ and z-i- Derivative of <fi blows up near 
these points. This is why we introduce scaling factor s. We can write a power 
series of <f> near z\ (near Z2 situation is the same by the symmetry) 



= c_ fe (z - Zl )- k +^ + c_ k+1 (z - Zl y k+ ^ + ■ • • + Co + . . . . 
This means that for s > 1 derivative can be estimated by 

| C _ fc |( S -l)- fc + 1 /2 + | C)b+1 | (s _ 1) - fc+ 3/ 2 + ... , 

Tail of this series can be estimated because series converges in a disc of a fixed 
radius (radius is \z% + 1|), and sum of tail can be estimated by a sum of a 
geometric progression. Writing these power series explicitly we find (for s — 
1.002) 

\(f>'\<55, \4>"\ < 11800, |</> (3) | < 8.69 x 10 6 , 

|^ (4) | < 1.08 x 10 10 , |^ (5) | < 1.90 x 10 13 , |</> (6) | < 4.25 x 10 16 , 
1^1 < 1.17 x 10 20 . 

The maximal values for first six derivatives of v are 

|i/|<0.28, K|<0-45, |z/ 3) |<1.12, 



(4) | < 3.69, |i/ 5) | < 15.3, |z/ 6) |<76.2. 
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The derivative dg\4>\ can be estimated by r\<j>'\. We can write sixth derivative 
of i/(|(/>|)|0'|/|0| as a rational function of partial derivatives of \<j)\, \<j)'\, and v. 
Than we apply triangle inequality and plug in the above estimates. Finally we 
have 



d9 e 



< 1.65 x 10 21 . 



Plugging e = 7r/5000 and estimate on sixth derivative into (fl"4"]) we find that 
error term in this case is less than 0.0034. 

Next we have to estimate modulus of continuity with respect to r. First we 
calculate 

d r \z - (a + bi)\ 2 = 2r- 2(acos0 + 6sin0). 



Applying this formula several times we find 

Iz-ll 



\<P\ 



Or 



r — cos o 
\z-l\S 



Ty/\Z - Zl\\z - 22 1 

r — x cos 6 — y sin 9 



< d r 



z-1 



2\z- Zl \ 2 S 



1 



s/\z - zi\\z - z 2 \ ^ 
r — x cos 9 + y sin 9 



2\z~z 2 \ 2 S 



z-1 



where S = y/\z — z\]}z — Z2\. Factoring out 

1 



2|z-l| ■ \z - Zl \ 5 / 2 \z - z 2 \ 5 / 2 



we get 



2(r - cos6»)|z - zx\ 2 \z - z 2 \ 2 - (r -xcosO - ysa\9)\z - l\ 2 \z - z 2 \ 2 

— (r — xcos9 + ys\a9)\z — l| 2 |z — zi| 2 
= -2(r 2 - l)(2cos 2 9r(x - 1) + cos 9{r 2 + l)(x - 1) + 2ry 2 ). 

This is a quadratic function with respect to cos#. Taking values of x and y into 
account we can write it as 

2n nf 1\ 1 592 

cos 9 + cos 9 \ r - 



r J 2 5625 

This quadratic function has two real roots. Their average is— (r + l/r)/2 < — 1, 
hence one root is definitely less than —1. The product of roots is a small negative 
number, which meant that the second root is positive and less than 1. Simple 
calculation shows that this root decreases as r grows. This means that the 
corresponding value of 9 increases. Hence it attains its maximal value at r = 1.4 
and the maximal value is at most 1.48. This gives us that the radial derivative 
of |</»'|/|(^| can be positive only on the arc 9 E [—1.48, 1.48]. By subharmonicity 
it attains the maximal on the boundary of {z | 1 < r < 1.4, —1.48 < 9 < 1.48}. 
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It is not very difficult to check that maximum is at z = 1.4 and it is equal to 
0.36. 
Let 



I(r) = / v(\Hr 1/k e l9 )\ 



<t>(r l l h e ie ) 



2tt' 



The derivative is 



^^(/;/(i^,Mf + /> (l ^(ffl)i 

By the symmetry the first integral is zero. In the second integral 



can be positive only when £ [—1.48, 1.48] and even in this case it is bounded 
by 0.36/(r 1 ^ 1 / fc fc27r). Hence 



I'(r) < 2 • 1.48 • 0.36/(r 1 - 1/fc /fc27r) < 0.0131 



If we compute values I(r{) and /(r^) (with precision 0.0034) then the minimum 
of P(y)/v on [r\, ra] is at least 



-J /fc (min{J(ri), I(r 2 )} ~ 0.0034 - 0.0131(r 2 - n) 



,i/fe-i 



(15) 



We take 3000 equidistributed points on [1, R] and compute I(r) at these points. 




20 30 40 50 60 
Figure 7: Values of I(r) 

The data for I(r) is shown on the Figure [7] Applying the error estimate fj 1 5[) 
we find a rigorous estimate from below of Pvjv. The minimum of Pv jv is at 
least 1.8079 which means that 

(3(1) > 0.2308. 
The Figure [8] shows the plot of log(fV/^)/logfc. 
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10 20 30 40 50 60 70 
Figure 8: Plot of \og{Pv / v) / logk 

5.4 Estimates of spectrum for other values of t 

We also computed lower bound on spectrum of the same snowflake for other 
values of t. Below are given base 13 logarithms of eigenvalues of discretized 
operator P (N = 1000, M = 400), lower bounds of \og(Pv/v)/ log k, t 2 /4 and 
upper bound of the universal spectrum from [5J [5] . For values of t close to zero 
we can not find a function that gives us a positive lower bound. 



t 


log 13 A 


m > 


t*/4 
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